% %south
dist_plus  = sum(reshape(f_plus_s,N_z,N_b),1);
dist_ss    = sum(reshape(f_dist_ss_s,N_z,N_b),1);
dist_0     = sum(reshape(f_0_s,N_z,N_b),1);
figure(44)
plot(log(supp_vec_s(1:N_z:end,1,1)),cumsum(dist_ss),'-b','LineWidth',6),hold on
plot(log(supp_vec_s(1:N_z:end,1,1)),cumsum(dist_0),'--r','LineWidth',6),
plot(log(supp_vec_s(1:N_z:end,1,1)),cumsum(dist_plus),':g','LineWidth',6),hold off
xlabel('$\ln b$','Fontsize',18,'Interpreter','Latex')
ylabel('Mass of firms','Fontsize',18,'Interpreter','Latex')
le=legend('Initial steady state','After exit','After exit and entry');
set(le,'Interpreter','Latex','Fontsize',16,'Location','NorthWest');
settt
grid on
% str1="figure_dist_0_s.pdf";
% str_save1=append(root,str1);
% settt
% print(gcf,'-dpdf',str_save1)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%north
dist_ss   = sum(reshape(f_dist_ss_n,N_z,N_b),1);
dist_plus = sum(reshape(f_plus_n,N_z,N_b),1);
dist_0    = sum(reshape(f_0_n,N_z,N_b),1);
figure(45)
plot(log(supp_vec_n(1:N_z:end,1,1)),cumsum(dist_ss),'-b','LineWidth',6),hold on
plot(log(supp_vec_n(1:N_z:end,1,1)),cumsum(dist_0),'--r','LineWidth',6),
plot(log(supp_vec_n(1:N_z:end,1,1)),cumsum(dist_plus),':g','LineWidth',6),hold off
xlabel('$\ln b$','Fontsize',18,'Interpreter','Latex')
ylabel('Mass of firms','Fontsize',18,'Interpreter','Latex')
le=legend('Initial steady state','After exit','After exit and entry');
set(le,'Interpreter','Latex','Fontsize',16,'Location','NorthWest');
%settt
% grid on
% str1="figure_dist_0_n.pdf";
% str_save1=append(root,str1);
% settt
% print(gcf,'-dpdf',str_save1)
